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Abstract. We study the spectrum of the Fibonacci Hamiltonian and prove 
upper and lower bounds for its fractal dimension in the large coupling regime. 
These bounds show that as A — > oo, dim{a-{H x)) ■ log A converges to an explicit 
constant (sj 0.88137). We also discuss consequences of these results for the 
rate of propagation of a wavepacket that evolves according to Schrodinger 
dynamics generated by the Fibonacci Hamiltonian. 



1. Introduction 

The Fibonacci Hamiltonian is a discrete one-dimensional Schrodinger operator 
[Hu]{n) = u{n + 1) + u{n - 1) + V{n)u{n) 
in The potential V : Z ^ R is given by 

(1) Vin) = \X[i-^-\i){n<j)-'+e mod 1), 
where A > is the coupling constant, is the golden mean, 

(2) , = + 



1 

1 + 

1 

1 



1 + 

and 9 E [0, 1) is the phase. 

This operator is important for both physical and mathematical reasons. On the 
one hand, it is the most popular quantum model of a one-dimensional quasicrystal, 
that is, a structure that shares many features with one displaying global order, 
but which in fact lacks global translation invariance. On the other hand, this 
operator has zero-measure Cantor spectrum and all spectral measures are purely 
singular continuous. These properties had been regarded as "exotic" in the context 
of general Schrodinger operators up until the 1980s, but for this operator family, 
they occur persistently for all parameter values. The Fibonacci Hamiltonian has 
been heavily studied since the early 1980s; see [7] for a recent review of the results 
obtained for it and related models. 

Let us recall some specific results and references that will be important for what 
follows. It is known that the spectrum of the Fibonacci Hamiltonian is independent 
of 9 (see, e.g., [4]). This follows quickly from strong convergence once one realizes 
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that for each pair 9, 9, there is a sequence rifc ^ oo such that + nk(p^^ converges 
to 9 in M/Z "from the right." We denote this common spectrum by S;^, 

Ex = (r{H) for every 61 e [0,1). 

It is natural to study the spectrum as a set. As wc shall sec, such a study is 
also motivated by the consequences one can draw for the long-time behavior of the 
solution of the time-dependent Schrodinger equation. 

It has been shown by Siito that the spectrum always has zero Lebesgue measure 

[35], 

(3) Leb(EA) = for every A > 0. 

This immediately implies the absence of absolutely continuous spectrum for all 
parameter values. ""^ It was later seen that one also has the absence of point spectrum 
for all parameter values; see Siito [34], Hof-Knill-Simon [17], and Kaminaga [21] for 
partial results and Damanik-Lenz [9] for the full result. Thus, the Fibonacci model 
exhibits purely singular continuous spectrum that is very rigid in the sense that it 
is not affected by a change of the defining parameters. 

The result on zero measure spectrum, (3), naturally leads one to ask about 
the dimension of this set. There are several popular ways to measure the fractal 
dimension of a nowhere dense subset of the real line. 

Let us recall the definition of two of these dimensions. Suppose we are given 
a bounded and infinite set 5 C M. A i5-cover of 5 is a countable union of real 
intervals, {/m}m>i, such that each of these intervals has length bounded by 5 > 0. 
For a € [0,1], let 

h"{S) = lim inf J2 

— >0 o-covcrs ^ — ^ 
m> 1 

It is clear that the limit exists in [0,cxd]. Moreover, if h°'{S) = for some a, then 
h" (S) = for every a' > a. Similarly, if h"{S) = oo for some a, then /i" (S) = oo 
for every a' < a. Thus, the following quantity is well-defined: 

dimH(S) = mi{a : h"{S) < oo} = sup{a : h"{S) = oo}. 

The number dim^j (S*) G [0, 1] is called the Hausdorff dimension of the set S. 

A different way to measure the fractal dimension of S is via the box counting 
dimension. The lower box counting dimension is defined as follows: 

dim^(5)=liminfi^i^, 

^ £^0 logl/£ ' 

where 

Ns{e) = #{i e Z : \js, (j + l)s) 5^0}. 
The upper box counting dimension, dim^(5'), is defined similarly, with the lim inf 
replaced by a lim sup. When dim^(S') and dim^(S') are equal, we denote their 
common value by dims (5) and call this number the box counting dimension of S. 
These dimensions are related by the inequalities 

dimH(5) < dimsC^) < dim+(5). 

In ReiK^ral, both inequalities may be strict; see, for example, [27, pp. 76-77]. 



Historically, these two properties were established in the reverse order. The methods used by 
Kotani [24] in the proof of absence of absolutely continuous spectrum (for almost every 0) were 
the key to proving zero measure spectrum. 
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The main goal of this paper is to study the fractal dimension of the spectrum 
of the Fibonacci Hamiltonian. The following result shows that for sufficiently large 
coupling, the dimensions just introduced coincide. 

Theorem 1. Suppose that A > 16. Then the box counting dimension ofT,\ exists 
and obeys 

dims (S a) = dimH(SA). 

While this theorem has not appeared in print explicitly before, it does follow 
quickly from a combination of known results. We present the relevant facts in the 
appendix. 

Theorem 1 is useful because it will allow us to obtain precise asymptotics for the 
fractal dimension of Sa as A ^ oo. The reason for this is the following. Tlw^ box 
counting dimension is easier to bound from below, while the Hausdorff dimension is 
easier to bound from above. Consequently, we will prove a lower bound for the box 
counting dimension in Section 3 and an upper bound for the Hausdorff dimension 
in Section 4. 

It is known how to describe the spectrum of the Fibonacci Hamiltonian in 
terms of the spectra of canonical periodic approximants. We will recall this in 
Section 2. The general theory of periodic discrete one-dimensional Schrodinger op- 
erators shows that the spectrum of a such a periodic operator is always given by a 
finite union of compact intervals. Our crucial new insight is a way to describe the 
asymptotic distribution of bandwidths in these periodic spectra. In this description, 
the following function plays an important role. Define 

f{x) = - [ (2 - 3a;) log 2 + (1 - a;) log(l - x) 

- (2a; - 1) log(2a; - 1) - (2 - 3a;) log(2 - 3a;)] 

on the interval (^,|)- Setting /(i) = log 2 and /(|) = 0, it is not hard to see 
that / extends to a continuous function on [5, |], and with the aid of symbolic 
computation one can confirm that it takes its maximum at the unique point 

12 - 2\/2 

X* = — ^ = 0.5395042867796 . . . , 

17 

with 

/* = f{x*) = log(l + ^/2)= 0.8813735870195 .... 



Write 
and 



Su{X) = 2A + 22 



Si{\) = i ((A - 4) + v/(A-4)2-12) 



With these functions of (sufficiently large) A we can now state the bounds on the 
fractal dimension of S;^ that we will prove in Sections 3 and 4, respectively. 

Theorem 2. (a) Suppose A > 4. Then 

r 



log5„(A)- 

(b) Suppose A > 8. Then 

r 

log 5; (A)- 



dimB(SA) > 



4 
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Since both Su{X) and 5';(A) behave asymptoticaUy hke log A, we obtain the fol- 
lowing result as an immediate consequence. We write dim for either dimjj or dims, 
which is justified by Theorem 1. 



Corollary 1. We have 



lim dim(I]A) •logA = /*. 

A— ^oo 



In particular, wc sec that the constant /* is the best possible in both bounds in 
Theorem 2. Let us compare our results with previously known ones. To facilitate 
this, we introduce 

f* = J— = ^°g(^ + ^) = 1.8315709239073 . . . , 
log (j) log (p 

so that our results can be summarized as follows. We have for A > 16, 
(4) f*T^. < di-B(S.) = dimH(E.) < /#- ^"S^ 



log5„(A) - ^' -V \ogSi{\y 

and therefore the asymptotic behavior is 

As wc will sec below, there are two competing scaling processes, one scaling with 
<j) (the Fibonacci numbers) and one scaling with A (the inverse of the width of a 
band in the approximating periodic spectra). Thus, it is natural to write a bound 
in the form "constant times l2i^" and then to optimize the constant. 

log A 

Raymond [33] proved an upper bound for dimjj (S;)^) that has a 2 in place of our 
/# in (4). A simplified version of our approach (which we will comment on later in 
the paper) quickly gives a lower bound with replaced in (4) by 1.5 and an upper 
bound with replaced in (4) by 2; the latter being Raymond's result. These 
numbers appear naturally in this context and are associated with the support of a 
certain discrete probability distribution. A more detailed study of this distribution 
then led us to the discovery of f^, which describes the actual asymptotic behavior 
of the fractal dimension of the spectrum as we saw above. 

Lower bounds for the dimension of the spectrum were initially obtained as a 
consequence of certain continuity properties of the spectral measures with respect 
to certain Hausdorff measures. Since the spectral measures are supported on the 
spectrum, one can obtain a lower bound for the Hausdorff dimension of the spectrum 
in this way. As mentioned above, this also bounds the box counting dimension from 
below by general principles. Spectral Hausdorff continuity results for the Fibonacci 
Hamiltonian were shown in [6, 8, 19, 22]. The best lower bound that has been 
obtained in this way can be found in [22] and it reads 

^'"^"^^^^ - i^Taxy 

where^ 

log (4^) 

K = / « 0.0126 

51og0 



■^Notice that there is a typo in [22]. They have k = log 2o''iog0 ' ^ negative number! 
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and 



C(A) = ^(logA + 0(l)). 

\0g(j) 



Thus, for A large, this gives a lower bound for dimiif(Ex) as in (4), but with 
replaced by 



(5) « 0.00188. 

6 log V 5 

Liu and Wen [25] then extended the approach employed by Raymond. They 
study the case of general frequencies. Specialized to the Fibonacci case, their result 
shows that for A > 20, 

log 2 . s log 3 



101og2 + 31og(4(A-8)) ^ ^^'"-(^^^ - log(A^- 
Let us discuss this result in the large coupling limit. Since 

log 3 

—2- « 2.28301 > 2, 
log0 

the upper bound does not improve Raymond's result. The lower bound has a 
constant coefficient 

log 2 

—2— « 0.48013 
31og0 

in front of a significant improvement over the result that can be extracted 

from [22]. Again, by our result, the optimal constant is /* « 1.83157. 

Our interest in obtaining the optimal constant docs not only stem from 
natural curiosity. An interesting and mathematically challenging problem is to 
study the spreading of a wavepacket in a quantum system in the case where the 
initial state has a purely singular continuous spectral measure. This is the case for 
every initial state from i"^ (Z) for a system governed by the Fibonacci Hamiltonian. 
One is often especially interested in the spreading of a wavepacket that is initially 
localized on just one site. 

That is, with H as above, we consider tp{t) = e~^*^6i and study its spreading 
via the time-averaged outside probabilities 

Pr.{N,T)=J2i e-^|(e-*^(5i,(5„)f 

n>N-^ Jo 

and 

Pi{N,T)= 5^ ^ re-^\{e-'*"S^,6r^)fdt. 

Let P{N, T) = Pi{N, T) + Pr{N, T) and define 

logP(r" - 2,T) 



S (a) = — lim inf ■ 

T — ^oo 

S'^{a) = — lim sup ■ 



T^cx) log T 

and 

log P(T" - 2,T) 



6 
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For every a, < S^{a) < S (a) < oo. These numbers control the power decaying 
tails of the wavepacket. In particular, the following critical exponents are of interest: 

= sup{a > : S^{a) < oo}. 

One can interpret as the rates of propagation of the fastest (polynomially small) 
part of the wavepacket; compare [15]. In particular, if a > q;+, then P(T",T) goes 
to faster than any inverse power of T, and if a > a~ , then there is a sequence of 
times Tfe ^ oo such that P(T^, T}..) goes to faster than any inverse power of T^. 

In Section 5 we prove a result for general Schrodinger operators on £^(Z) that 
will imply the following consequence for the Fibonacci Hamiltonian. 

Theorem 3. For every A > and every 6 G [0, 1), we have that 

at > dim|(SA). 

Consequently, for A > 4 and every 6, we have 



To discuss this result in the large coupling limit, let us be slightly inaccurate^ 
and write 

There are two main previous approaches to quantum dynamical lower bounds for 
the Fibonacci model. The first is based on spectral continuity and the papers 
[6, 8, 19, 22] contain results obtained in this way. For = 0, the best bound 
is contained in [22] and it has (6) with « 1.83157 replaced by (5), that is, 
w 0.00188. For other values of 6, the best bound can be found in [8] and the 
constant in this boimd is even smaller. 

The other approach is based on complex energy methods and the Plancherel 
Theorem; see [11] (and also [12] for a way to combine the two approaches). For 
6 = 0, the paper [11] has (6) with replaced by 1/6. It is possible to treat general 
6 along the same lines using [10], but the dynamical lower bound has a somewhat 
smaller constant in the general case. 

Thus, on the one hand, our result improves the constant from the previously 
best value 1/6 to ~ 1.83157 and, on the other hand, this is the best one can 
do using the method put forth in this paper. We would like to mention that [13] 
contains the following upper bound for a+, which holds for A > 8, 

at < 2.^^ 



While we know that the dimension of the spectrum indeed behave like /# in 
the large coupling limit, we expect that the dynamical quantities at behaves like 



log Si{\) 

;he spectri 
t the dyna 

2 1^1^ in the large coupling limit. That is, we expect the following to hold, 

hm a„ 1 = 2. 

A— ►oo log (p 

The reason for this is that the spreading of the fastest part of the wavepacket is 
determined by the "most continuous" part of the spectral measure and, in this case. 



^The precise statement is that liminfx_,oo ctu IHg^ — 
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by the region in the spectrum that is "the thickest." Since we will sec that there 
is indeed a small region that is thickest in a natural sense, the factor 2 will then 
appear naturally. A forthcoming publication, [14], will deal with this issue using 
ideas and results from [13]. 

2. The Band-Gap Structure of the Approximating Periodic Spectra 

In this section wc describe the canonical coverings of T,\ by (unions of) periodic 
spectra and the hierarchical structure of these sets. We use the combinatorics of 
this description to derive detailed results about the distribution of bandwidths in 
these spectra. 

For E & R and A > 0, we define a sequence of numbers Xk = Xk{E, A) as follows. 

(7) x-i=2, xq=E, xi = E-X, Xk+i = XkXk-i - Xk-2 for fc > 1. 

Using this recurrence and the initial values, one can quickly check (see Siito [34]) 
that 

(8) ^k+i + ~ Xk+iXkXk-1 = 4 + A^ for every fc > 0. 
For fixed A > 0, define 

ak = {E€^:\xk{E,X)\ < 2}. 

The set Uk is actually equal to the spectrum of the Schrodinger operator H whose 
potential Vk results from V (with ^ = 0) by replacing with Fk-i/Fk (cf. [34]). 
Here, {Fk}k>o denotes the sequence of Fibonacci numbers, that is, 

Fo = 1, Fi = 1, Ffe+i =Fk+ Fk-i for fc > 1. 

Hence, Vk is periodic and Uk consists of Fk bands (closed intervals). Siito also 
proved that 

(9) (cTfc-l U CTfc) D (cTfc U (Jk+l) 

and 

(10) = fl ((Tfc U Cfc+i) = {E : {xk} is a bounded sequence}. 

fe>i 

Prom now on, we assume 

(11) A>4, 

since we will make critical use of the fact that in this case, it follows from the 
invariant (8) that three consecutive XkS cannot all be bounded in absolute value 
by 2: 

(12) cTfc n CTfc+i n (Jk+l = 0. 

The identity (12) is the basis for work done by Raymond [33]; see also [13, 22], 
which describe the band structure on the various levels in an inductive way. Let 
us recall this result. Following [22], we call a band Ik C ak a. "type A band" if 
Ik C <Jk-i (and hence Ik O {(Tk+i U (Tk-2) = 0)- Wc call a band 7^ C (T/c a "type B 
band" if Ik C ak-2 (and therefore Ik n ak-i = 0). Then, we have the following 
result (Lemma 5.3 of [22], essentially Lemma 6.1 of [33]): 
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Lemma 1. For every A > 4 and every k > 1, 

(a) Every type A band Ik C ak contains exactly one type B hand Ik+2 C (Tk+2, and 
no other bands from ak+i, (7k+2- 

(b) Every type B hand Ik C Uk contains exactly one type A hand Ik+i C Ck+i and 
two type B bands from ak+2, positioned around Ik+i- 

Lemma 2. For every band I of Uk, we have that I n Y,\ ^ 0. 

Proof. Let / be a band of ak- Choose a band /^^^ in Gk+i U ak+2 with / D I^^\ as 
is possible by Lemma 1. Iterating this procedure, we obtain a nested sequence of 
intervals and hence a point E G I for which the corresponding trace map orbit is 
bounded. Thus, by (10), E €i:x. □ 

We define 

ak = number of bands of type A in o-fe, 
bk = number of bands of type B in o-fc. 

By Raymond's work, it follows immediately that ak + bk = Fk for every k. In 
fact, we have the following result. 

Lemma 3. The constants {ak} and {bk} obey the relations 

(13) ak = bk-i, bk = ak-2 + 2bk-2 
with initial values 

ao = 1, ai = 0, bo = 0, bi = 1. 
Consequently, for k>2, 

(14) ak = bk-i = Fk-2- 

Proof. The recursions (13) hold by definition. The explicit expressions in (14) then 
follow quickly by induction. □ 

We axe interested in the size of the bands of a given type on a given level. To 
capture the distribution of these lengths, we define 

ak,m = number of bands b of type A in Uk with #{0 < j < k : bdaj 7^ 0} = m, 
bk,m = number of bands b of type B in ak with #{0 < j < k : b (1 aj ^ ^} = m. 

The motivation for this definition is the following: If we consider a given band and 
its location relative to the bands on the previous levels, each time the given band 
is contained in a band on a previous level, we essentially pick up a factor roughly 
of size A"'^ (for A large). Thus, for example, there are ak.m bands of size w A"™ 
of type A in ak- Since the combinatorics of the situation is A- independent, we 
choose to separate the two aspects. This will allow us to get much more precise 
information in the regime of large A, which is extremely hard to study numerically. 

Lemma 4. We have 



(15) 



ak,m = bk-l,m-l, bk,m = afe-2,m-l + 26fe_2,m-l 
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with initial values 



(io,m = for m > and ao,o = 1, 

ai,m = for m > 0, 

bo,m = for m > 0, 

bi,m = for m > and bifi = 1. 



Consequently, 

(16) ak,m — bk-l,m-l 



otherwise. 



Remark. The fact that Ofe.m is zero when [|] < m < [^J fails is an immediate 
consequence of the properties (9), (10), and (12) established by Siito. As a slight 
variation of the proof of Theorem 4 below shows, this fact alone is sufBcient to give 
a quick proof of 



Our more detailed description of the numbers ak,m will then enable us to prove the 
stronger lower bound with /# in place of 1.5, which is optimal as discussed in the 

Introduction. 

Proof. It is immediate from the definition and (13) that the recursions (15) hold. 
Recall that the Chebyshev polynomials of the first kind are defined by the recurrence 
relation 

To(x) = 1, 
Ti{x) = X, 

Tm+l{x) = 2xTm{x) - Tm-l{x). 

Write Tm{x) as 

^m(*^) ~ ^ ^ ( 1) ^r,mX 

On the one hand, it is known (see, e.g., [1]) that for m > 1 and < r < [^Ji 

(17) Cr,m = 2^-'^-'^^h~'' 

m — r \ r 

On the other hand, the recursion generating the polynomials says that Co,o = 1) 
Co,i = 1, and 



^m—2r—l 



£ (-l)'^c,,„+ix'"-2'-+i = 2x£(-l)'-c,,„x™-2'-- g (-l)'-c,,„-i 

r=0 r=0 r=0 

= g(-l)'^2c,,„a;— 2'-+!+ (-l)'^+'c,,„_ia; 

r=0 r=0 

= ^(-l)'^2c.,„a;'"-2'-+i+ {-iycr-i,m-ix"'-^''+'. 

r=0 r=l 



m—2r—l 
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It follows that 

(18) Cr,m+1 — '^Cr,m ~t~ Cr—l,m—l' 

Denote dk,m = C2m-k,m- Then, using (18), we find that 

afe+l,m+l = C2(TO+l)-(fe+l),m+l 
= C2m-fe+l,m+l 
— 2c2to— fe+l,m "I" C2to— fe,m— 1 
= 20fc_i,m + a-k-2,m-l- 

Comparing this with the recursion 

flfc+l,m+l — '2ak-l,m + a*;-2,m-l 

we established above, along with the initial values, we see that 

(19) Clk^m — dk,m — ^2m—k,m- 

Combining (17) and (19), we obtain (16). □ 
Recall that we introduced the function 

/(x) = - [ (2 - 3a:) log 2 + {1 - x) log(l - x) 

- {2x - 1) log(2x - 1) - (2 - 3a;) log(2 - 3a;)] 
on the interval (i, |). We set /(^) = log2 and /(|) = and write 

X* = — = 0.5395042867796 . . . 

17 

for the unique point in (|, |) where / takes its maximum value, 

/* = f{x*) = log(l + \/2) = 0.8813735870195 .... 
Proposition 1. For ^ < m < we have 

(20) fc- V2 exp (m/ (^)) < ak,m < k''^ exp (m/ (^)) . 
Consequently, 

(21) lim max — loga^t = /*. 

fc— ►oo m rn 

Proof. As we saw above, when | < m < 

„ _ 92fc-3m-i f k-m 



k — m \ 2m — k 



Thus, we see that a^ k = 22 ^ and a^. 2k = 1, so the estimate (20) holds when 



m = I or m = ^. 

Let us now consider | < m < ^ , in which case 

(22) -<k-m<-, 

k 

(23) l<2m-k<-, 

(24) l<2fc-3m<|. 
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Stirling's approximation gives 

for every n > 1, where we write a>^b if a <b and a>b. Thus, 

k-m (fe^)^— 



2fe-3m 



k - m \{2m - k){2k - 3m) J ^ 2m-fc ^2"'-fe ^ 2fc-3m ^ 
Due to the estimates (22)-(24), (20) will follow once we show that 

Writing x = m/k^ we find that 

/7 \k—m 

^ (2m -fc)" (2fc-3mf'=-^'" 

_ 22fc-3fex (fc - fca;)^ 

{2kx - (2fc - 3A;ar)''=-^'=" 



2k— Sm ' 



= exp {kxf (x)) 



= exp(m/(^)). 



This finishes the proof of (20), which in turn immediately implies 

lim sup max — log am k < /* • 

On the other hand, as k gets large, we can choose m so that m/k gets arbitrarily 
close to X*, so that by (20) again, 

lim inf max — log Um k > /* • 

fe— >C!0 m m ' 

Thus, we have established (21). □ 

3. A Lower Bound for the Box Counting Dimension of the Spectrum 

In this section we will prove a lower bound for the (lower) box counting dimension 
of the spectrum of the Fibonacci Hamiltonian. Recall that the lower box counting 
dimension of a bounded set 5 C M is defined as follows: 

dimo(5) = lim inf ^sis) 
J logl/£ ' 

where 

Ns{e) = #{i e Z : \je, {j + l)e) n 5 ^ 0}. 
The following was shown in [11]. 
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Lemma 5. For every A > 4, there exists Su{X) such that the following holds: 

(a) Given any {type A) band Ik+i C ak+i lying in the band Ik C ak, we have for 

every E eh+i, 



< Su{\). 



(b) Given any {type B) band Ik+2 C (7k+2 lying in the band Ik C ak, we have for 
every E e Ik+2, 

^k+2iE) 



< 5„(A). 



For example, one can choose S'„(A) = 2A + 22. 

We use the symbol Su (A) instead of the explicit 2A + 22 to make the dependence 
of everything that follows on this quantity explicit. 

Inductively, Lemma 5 gives an upper bound for on each band of a/, and this 
in turn will give a lower bound on the length of a band, which will be crucial in the 
proof of the following result. 

Theorem 4. For all A > 4, 

(25) dimB(SA) > ^* 



log5„(A) 

Proof. Let mk = [3kx*\. Since x* > 0.4, the sequence {ruk} is strictly increasing 
and limfe^oo ^ = x* ■ Write 

fk = — log a3fe,mfe • 

mfe 

We know from Proposition 1 that 

(26) hm fk - /* 

and it is obvious that 

(27) hm = 1. 

fe— >(x) mfe 

For a given k. let us consider the type A bands in (73fe that lie in bands on 
previous levels. By definition, there are Nk := a^k,mk such bands. By Lemma 5, 
each of them has length at least efc := 45„ (A) since on each band of crsfe, the 
function x^k is strictly monotone and runs from ±2 to =p2. Let {^3fej}^'i be these 
bands, enumerated so that A^k,] is to the left of A^k^j+i for every j. Each band 
has non-empty intersection with by Lemma 2. Thus, for every j, there exists 
ETik.j G ^3fej n Ea- Clearly, the energies {E^k.j} are increasing in j. Consider 
{Ezk,j} with j odd, that is, j = 2s + 1, < s < [_Nk/2\. Since any two bands 
^3A:,2s-i, ^3fc,2«+i are separated by the band A3k,2s, which has length at least efe, 
we get \x^k.2s-i - a;3fc,2s+i| > £fe for every s. Thus, the E^k.2s+i belong to different 
£-boxes if £ < £fe. We conclude that N-£^{e) > ^ for every £ < £fe. 
Given any £ > 0, choose k with £fc+i < £ < £fc. Then, 

log A^E, (£) ^ log(a3fc,™, ) - log 2 ^ fk- ^\og2 



> — ^^^^^ '—^ — > 



log l/£ - log l/£fe+i - ^ log 5„(A) • 
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Since, as e ^ 0, we have k, mk 



lim 



oo and 
Ik 



r 



fc^oo I^log5„(A) log5„(A)^ 
by (26) and (27), the result follows. 



□ 



4. An Upper Bound for the Hausdorff Dimension of the Spectrum 

In this section we prove an upper bound for the Hausdorff dimension of S^. We 
will use the canonical coverings dk U <Jk+i of T,\ and estimate the lengths of the 
intervals in those periodic spectra from above using our combinatorial result from 
Section 2 together with a scaling result that is analogous to Lemma 5, but which 
gives bounds from the other side. 

Namely, the following was shown by Killip, Kiselev, and Last [22, Lemma 5.5]. 

Lemma 6. For every A > 8, there exists Si{\) such that the following holds: 

(a) Given any {type A) hand Ik+i C ak+i lying in the hand Ik C ak, we have for 
every E € Ik+i, 

^k+iiE) 



^kiE) 



> Si{X). 



(b) Given any {type B) hand Ik+2 C <Jk+2 lying in the band Ik C ak, we have for 
every E e 4+2, 



{E) 



> Si{X). 



For example, one can choose 

Si{X) = 



^((A-4) + V(A-4)2-f2). 



As before, we use the symbol Si (A) instead of the explicit possible choice to make 
the dependence of everything that follows on this quantity explicit. 

Inductively, Lemma 6 gives a lower bound for \x'^.\ on each band of ak and hence 
an upper bound for the length of a band. 



Theorem 5. Suppose A > 8. Then 

dimH(SA) < 



log Si{X)- 



Remark. Raymond [33] proved an upper bound of the form 

21og(/) 



dimH(SA) < 



log Si{X)- 



As above, to get this weaker result, all one needs is information about the support of 

o-k.rrn which in turn follows quickly from the properties (9), (10), and (12). Using our 
more detailed information about the values of ak,m on the support, we can improve 
the constant 2 to the value « 1.83157, which is optimal by our discussion in 
the introduction. Our proof is inspired by Raymond's proof, and the improvement 
stems from our more detailed analysis of the distributions of bandwidths in the 
approximating periodic spectra. 
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Proof. As recalled above, the set (JkUak+i is a finite union of compact real intervals 
that covers Ea- There are ak,m + bk,m bands in at that lie in exactly m intervals on 
previous levels. We know that the length of each of these bands is bounded from 
aboveby 4S',(A)-'". 

Thus, it suffices to show that, given any 

we have 

LfJ L^J 
lim V {ak,m + bk,m){^Si{X))-'"'+ V (afc+i,„ + 6fc+i,„) (45KA))-«™ = 0. 

fc— >oo 

For simplicity, we will only consider 

Lf J 

the other terms can be dealt with in an analogous way. 

By our upper bound for ak,m established earlier and the assumption (28), we 
have 

Ck<k'/' E exp(m(/(f)-slog5;(A))) 

m=r|i 

<k'/'cxp Q(r-5log5KA)) 

Using (28) again, we see that the right-hand side goes to zero as fc ^ oo. □ 
Theorem 4 and Theorem 5 together establish Theorem 2 from the Introduction. 

5. Transfer Matrices, Box Counting Dimension, and Wavepacket 

Spreading 

In this section we will first consider general half-line Schrodinger operators and 
derive a number of consequences from polynomially bounded transfer matrices. We 
then turn to general Schrodinger operators on the line and use our half-line results 
to derive a lower bound for the exponent governing the spreading of the fastest part 
of the wavepacket in terms of the box counting dimension of the spectrum. These 
results are relevant in our context since for the Fibonacci potential, it is known 
that the transfer matrices are polynomially bounded for all parameter values and 
all energies in the spectrum. 

Consider a discrete half-line Schrodinger operator on £^(Z+), 

(29) [H+u]{n) = u{n + 1) + u{n - I) + V{n)u{n), 

with Dirichlet boundary condition u(0) = 0. Here, V is a bounded real-valued 
function. if+ is a bounded self-adjoint operator in €^(Z+). The vector 5\ is cyclic 
and we denote its spectral measure by ^. 

For z G C and m, n G Z+, we denote by T(n, m; z) the transfer matrix associated 
with the difference equation 

(30) u{n + u{n + V{n)u{n) = zu{n). 
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That is, T(n, m; z) is the unique unimodular 2x2 matrix for which we have 

u{n + l)\ , .fu{m+l) 
^ , . ' ] = T[n,m;z) / x 
u{n) J \ ' w y u{m) 

for every solution of (30). 

To study the spreading of a wavepacket with initial state Si under the dynamics 
generated by H, one usually considers for p > 0, 

One is interested in the power-law growth of this quantity and therefore studies the 
lower transport exponent 

Hs,\yj plogT 
and the upper transport exponent 

pI [p] = hmsup — 7- . 



Both functions {p) are nondecreasing in p and hence the following limits exist: 

Alternatively (see [15, Theorem 4.1]), we also have = sup{a > : 5'^ (a) < oo}, 
where 

5-(a) = -liminf^"g^(^"-^'^) 

^ ' T^oo logT 

.™„„„iog^+(r"-2,r) 



S'^ (a) = — lim sup 



and 



9 f°° 

^ Jo 



n>N 

The following result derives several consequences from the fact that the transfer 

matrices arc polynomially bounded for energies from some subset of the real line. 
We denote by /i the spectral measure associated with the vector 5i. 

Theorem 6. Consider a half-line Schrddinger operator with bounded potential 
as above. Let A c [—.8,-8], B > 0, and assume that there exist positive constants 
C, a such that for every E G A and every N >1, 

(31) \\T{n,m;E)\\ < CN" forallm,n& [1,7V]. 

(a) For every cr > 0, there exists £0(0") > such that 

(32) n{[E-e,E + e\)>De^ 

for e € (0, So) o,nd E £ A. Here, L = {1 + a){l + 3a) and D depends only on a. 

(b) For any p > 0, we have that 

(33) - ' + 
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In particular, 

(34) at > dim|(A). 

Proof. It was shown by Germinet, Kiselev and Tcheremchantsev [15, Proposition 
2.1] that for any M > 0, 

j-E+e/2 

(35) ^i{[E-e,E + e])>Cl \\T{N,Q-x)\\-'^ dx - €28^ . 

J E-e/2 

Here e e (0,1), E e [-B,B], N = [e~'^~''\, cr > 0, and the constants Ci,C2 
depend only on B,M,a. 

Let E & A. To bound the integral in (35) from below, we need an upper bound 
for ||T(A'', 1; x)|| for x close to E. This can be accomphshed using a perturbative 
argument of Simon; compare [11, Lemma 2.1]. Namely, if 

D{N,E)= sup \\Tin,m;E)\\, 

l<n,m<N 

then for any (5 G C and 1 < n < we have 

\\T{n,l;E + d)\\ < D{N, E) cxp{D{N, E)\n\S). 
By assumption, D{N, E) < CN" ioi E e A and N>1. Thus, 

||r(7V,0;i; + (5)|| < CiV"exp(C7Vi+"(5), 
and for any 5 with jjj < 7 = 1/2N~^~", we get 

(36) \\T{N,1;E + S)\\<C'N" 

with a uniform constant. 

Since N = le~^~'^\, e G (0,1), wc see that 7 < e/2 for e small enough and 
therefore [E - E + j] C [E ~ e/2, E + e/2]. Using (35) and (36), wc find 

n{[E-e,E + e])>Ci / \\T{N,0;x)\\-^dx - C2e^ 

> C377V-^" - C2£^ 
= C47V-^-^" - C2e^ 

with appropriate constants. Taking M = (1 + (t)(1 + 3a) + 1, wc obtain (32). 

In order to prove the dynamical lower bound (33), we will apply [37, Corol- 
lary 4.1]; see also [2]. In our notation, this result reads 

(37) /3±(p)>7?±(g), (^,1), 

where D^{q), < g < 1 are the generalized fractal dimensions of the spectral 
measure; compare [2, 3]. One of the equivalent ways of defining these dimensions 
is the following (cf. [37, Theorem 4.3]): 

(38) !,+ (,) =linisup;^^i%£l, 

e^o ((?-l)loge 

and similarly for D~{q) with limsup replaced by liminf, where 
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To bound Sfj,{q, e) from below, wc can follow the proof of [37, Theorem 4.5]. Let 
Ij = [je, (j + Denote by Ja the set of j € Z for which Ij D A ^ 0. Thus, for 
every j € Ja, there exists Ej e [je, {j + l)e) such that Ej e A. The bound (32) 
impUes for e small enough, 

+ Klj) + > ^^{[E, - e, Ej + £]) > De^, 

for every j € Ja, where L = (1 + (t)(1 + Sa). Thus, 

On the other hand, since (a + 6)* < a* + 6^, a, 6 > 0, g e (0, 1), one has B < 
3Sf_i{q,e). Therefore, 

(39) S^{q, s) > C{q)e'^^ ^ 1 = C{q)s^'^N{e) 

with 

N{e) = #{j€Z : IjnA^H)}. 
It follows from (38)-(39) that 



(40) D^(.)>^^^, ..(0,1). 

The bounds (37) and (40) imply, for any q G {j^, 1), 

^ dnn^^jA) - qL _ 

Since this holds for L = {l + a){l + 3a) with any a > 0, letting a ^ and q j^, 
we obtain (33). 
Since 

aj= lim PZ(p), 

p— >+(X) ^ 

the estimate (34) follows. □ 

Remark. One can improve the bounds (32), (33) if one has a nontrivial (i.e., better 
than ballistic) upper bound on a+. In this case [38], (35) holds with N = [e^^'J , 
where p = + a, a > 0. Thus, (32) holds with a smaller value of L. 

Let us now turn to the whole- line case and prove a result analogous to (34). 
Consider the operator H acting on £'^{Z) as 

[Hu]{n) = u{n + 1) + u{n - 1) + V{n)u{n) 

with bounded F : Z ^ M, and choose K so that a{H) C [-K + l.ii' - 1]. The 
transfer matrices T(n, m; z) are defined in an analogous way. 

For N > 1, recall that the time-averaged right outside probabilities are given by 



P,(iV,r) = |j^°°e-*^|(e-*^5r,<5„)| 



2 , 

dt. 



n>N 

It was proved in [13] that 

(41) CT-^Ir{N, T) < Pr{N, T) < Ce-"^ + CT^Ir{N, T), 
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where 

Ir{N,T)= r \\T{N,l-E+^)\\-^dE. 

J-K 

Although not stated cxphcitly, this result from [13] can be extended easily to the 
half-line case (29). Thus, for the time-averaged outside probabilities P+{N,T), we 
have 

(42) CT-^I+{N, T) < P+{N, T) < Ce'"^ + CT^I+{N, T), 

where 



I+{N,T)= [ \\T{N,l-E+^)\\-^dE. 

J-K 



Corollary 2. Suppose that H is a discrete Schrddinger operator on the line with 
hounded potential. Assume A<z[—B,B](Z'& is such that for E G A and N > 1, 

(43) \\T{n, m; ^)|| < CN" for all m,nG [1, N] 

with uniform constants C,a. Then, for the initial state di, we have that 

>dim±(^). 

Proof. Let if+ be the half-line operator of the form (29) with potential V+ given 
by V+(n) = V{n) for n > 1. The transfer matrices are the same for both operators 

since they are defined locally, and thus wc have Ir{N,T) — I+{N,T). The condition 
(31) holds for the operator and Theorem 6 therefore yields 

a±+ > dim±(A) 

for the corresponding dynamics. 

Let us assume that dimg(^) > (if dim^(A) = 0, the result is trivially true) 
and show that a~ > dim^(A). Take any a with < a < dim^(A) < a~_|_. Due to 
the definition of a^.,., we have that 



(44) P+(T" -2,T) > T" 



-M 



with some finite M > for T sufficiently large. Since Ir{N,T) = I+{N,T), it 
follows from (41) and (42) that 

PriN,T) > CT-^I+{N,T) > CT-^ {P+{N,T) - Ce'''^) . 

It follows from (44) that for T sufficiently large, 

Pr(r" - 2,T) > CT-^-^. 

Consequently, for the full time-averaged outside probabilities, we find 

P{T--2,T) = ± e-^ J2 \{e-''"S,,6n)\' dt 

^ -^0 |n|>T» 

> CT-^-^. 

It now follows directly from the definition of a~ that a~ > a. Since this holds 
for every a < dimjj(A), the result follows. For a+, the proof is the completely 
analogous. (The bound (44) then holds for some sequence of times.) □ 
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Remarks, (a) One can prove a similar result under the condition 
\\T{n,m;E)\\ < CN" for all m,n G [-N,-l], 
where E e A, N > 1. 

(b) This establishes Theorem 3 as stated in the Introduction since it has been shown 
that in the Fibonacci case, the transfer matrices are polynomially bounded for all 
energies in the spectrum with uniform constants C, a. See lochum-Testard [18] for 
the case 6 = and Damanik-Lenz [10] for the case of general 6. 

Appendix A. The Hyperbolicity of the Trace Map at Large Coupling 

AND Some of its Consequences 

A.l. Description of the Trace Map. The main tool that we are using here is 
the trace map. It was originally introduced in [20, 23]; see also [34] for proofs of the 
results described below. Let us recall that the numbers Xk = Xk {E, A) introduced 
in Section 2 satisfy the recursion relation 

(45) Xk+i = XkXk-i - Xk-2, 

with initial conditions X-i = 2, xo = E, x\ = E — X, and the invariance relation 

(46) xl_^_y - Xk+iXkXk-i = 4 + 

for every k G Z_|_. Because of (45), it is natural to consider the so-called trace map, 

T-.R^^R^, T{x, y, z) = {xy - z, x, y). 

The sequence {xi, 0:2, 2:3, . . .} can be considered as the sequence of first coordinates 
of points in the trace map orbit having initial condition (a;i, a;o, a;-i). 
By (46), the following function is invariant under the action of T: 

I = x^ + y"^ + z"^ — xyz — 4. 

In other words, T preserves the family of cubic surfaces 

S/ = {{x,y,z) e : x'^ + y'^ + z'^ - xyz - 4 = I}. 

The surface So is called the Cayley cubic. Denote by l\ the line 

lx = {{E-\E,2) -.EeR). 

It is easy to check that l\ C • 

The following result, proved in [34], characterizes the spectrum of the Fibonacci 

Hamiltonian in terms of trace map dynamics and therefore establishes an important 
and fruitful connection between spectral and dynamical issues in this context. 

Theorem 7. The energy E belongs to the spectrum, of the Fibonacci Hamilton- 
ian if and only if the positive semiorbit of the point {E — A, E, 2) under iterates of 
the trace map T is bounded. 
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A. 2. Hyperbolicity of the Trace Map for Large A. Denote by f\ the restric- 
tion of the trace map to the invariant surface S\2 . That is, f\ : S\2 8^2 and 
fx = ^|S;^2- Denote by the set of points in §;^2 whose full orbits under fx are 
bounded. 

Let us recall that an invariant set A of a difFeomorphism / : M ^ M is locally 
maximal if there exists a neighborhood U (A) such that 

A = f|, /"([/). 

An invariant closed set A of a diffeomorphism / : M — > M is hyperbolic if there 

exists a splitting of the tangent space T^M = ® EJ^ at every point x E A 
such that it is invariant under Df, and Df exponentially contracts vectors from 
the stable subspaces {E^} and exponentially expands vectors from the unstable 
subspaces {E^}. 

See [16] for a detailed survey of hyperbolic dynamics and an extensive list of 
references. 

Casadgli proved the following result [5] . 

Theorem 8. For every A > 16, the set Qx is a locally maximal invariant hyperbolic 
set of fx ■■ §A2 • 

A. 3. Some Properties of Locally Meiximal Hyperbolic Invciriant Sets of 
Surface DifFeomorphisms. Consider a locally maximal invariant transitive hy- 
perbolic set A C M, dimM = 2, of a diffeomorphism / G Diff''(M), r > 1. 
We have A = n„ez /"(C^(A)) for some neighborhood f/(A). Assume also that 
dim£;" = dimE' = 1. 

Let us gather several results in this general context that we will eventually spe- 
cialize to the case; where A is given by fix for fx '■ ^ ^A^, ^ > 16. 

I. Stability. There is a neig hborhood U C Diff^(M) of the map / such that for 
every g gU, the set Ag = n„gzff(t^(A)) is a locally maximal invariant hyperbolic 
set of g. Moreover, there is a homeomorphism /i : A — > Ag that conjugates /|a and 
g\Ag, that is, the following diagram commutes: 



A 


/Ia 


> A 


h 








h 












A. 


> Ag 



II. Invariant Manifolds. For a; e A and small s > 0, consider the local stable 
and unstable sets 

W^{x) = {weM: d(/"(a;), /"H) < £ for all n > 0}, 

W^(x) = {wGM : d{f"(x), /"(w)) < e for all n < 0}. 

If £ > is small enough, then these are embedded C""-disks with TxW^{x) = E^ 
and TxW^{x) = EJ^. Define the (global) stable and unstable sets as 

W^{x) = U„ez+/-"(W^,^(x)), W^{x) = U„ez+r(W«(a;)). 

Define also 



W'{A) = U^^aW%x) and W{A) = U^eA^ix). 
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III. Invariant Foliations. A stable foliation for A is a foliation of a neighbor- 
hood of A such that 

(a) for each x G A, T{x), the leaf containing a;, is tangent to i?^; 

(b) for each a; S A, sufficiently near A, f{!F^{x)) C !F'^{f{x)). 
An unstable foliation .F" can be defined in a similar way. 

For a locally maximal hyperbolic set A C M of a C^-diffcomorphism / : M M, 
dimM = 2, stable and unstable C° foliations with C^-leaves can be constructed 
[26]. In the case of C^-diffeomorphisms, invariant foliations exist (see [30], 
Theorem 8 in Appendix 1). 

IV. Local HausdorfF Dimension and Box Counting Dimension. Consider, 
for a; G A and small e > 0, the set W"{x) fl A. The Hausdorff dimension of this set 
does not depend on a; e A and e > 0, and coincides with its box counting dimension 
(see [28, 36]): 

dimHW^{x) n A = dims M^" (a:) n A. 

In a similar way, 

dimnW^ix) n A = dims W/ (a;) n A. 
Denote ft.^ = dimHW^{x) n A and h" = dimHW;^{x) A. We will call and /i" 
the local stable and unstable Hausdorff dimensions of A, respectively. 

V. Global Hausdorff Dimension. Moreover, the Hausdorff dimension of A is 
equal to its box counting dimension and 

dimjj A = dims A = + /i"; 

see [28, 31]. 

VI. Continuity of the Hausdorff Dimension. The local Hausdorff dimensions 
/i''(A) and /i"(A) depend continuously on / : M ^ M in the C^-topology; see [28, 
31]. Therefore, dimjj Aj^ = dim^ Aj = h^{Af) + /i"(Aj) also depends continuously 

on / in the C"'^ -topology. Moreover, for r > 2 and C""-diffeomorphisms f : M ^ M, 
the Hausdorff dimension of a hyperbolic set Ay^ is a C'^~^ function of /; sec [26]. 

Remark. For hyperbolic sets in dimension greater than two, most of these properties 
do not hold in general; see [32] for more details. 

A. 4. Implications for the Trace Map and the Spectrum. Due to Theorem 8, 

the properties I- VI can all be applied to the hyperbolic set fix of the trace map 
fx ■ S;^2 for every A > 16. One can extract the following statement from the 

material in [5, Section 2]. 

Lemma 7. For A > 16 and every x eflx, the stable manifold W''{x) intersects the 
line Ix transversally. 

The existence of a C-'^-foliation JF" allows us to locally consider the set {il.x)C)lx 
as a -image of a set W"{x)r\Qx- Therefore, we obtain the following consequences 
for the spectrum of the Fibonacci Hamiltonian. 

Theorem 9. For A > 16, the following statements hold: 

(a) The spectrum Ea depends continuously on A in the Hausdorff metric. 

(b) We have dimHi^x) = dims {T,x). 

(c) For every small e > and every E gT,x, we have 

dimn {{E - s, E + e) nT,x) = dimH{T,x) 
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and 

dims {{E -s,E + e)nEx) = diniB(SA). 
(d) The Hausdorff dimension dim.H{T,\) is a C°° -function of X. 

In particular, part (b) establishes Theorem 1 as formulated in the Introduction. 
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